A 3D Vector-Additive Iterative Solver for the Anisotropic Inhomogeneous Poisson Equation in the Forward EEG problem
نویسندگان
چکیده
We describe a novel 3D finite difference method for solving the anisotropic inhomogeneous Poisson equation based on a multi-component additive implicit method with a 13-point stencil. The serial performance is found to be comparable to the most efficient solvers from the family of preconditioned conjugate gradient (PCG) algorithms. The proposed multicomponent additive algorithm is unconditionally stable in 3D and amenable for transparent domain decomposition parallelization up to one eighth of the total grid points in the initial computational domain. Some validation and numerical examples are given.
منابع مشابه
A 3D Finite-Difference BiCG Iterative Solver with the Fourier-Jacobi Preconditioner for the Anisotropic EIT/EEG Forward Problem
The Electrical Impedance Tomography (EIT) and electroencephalography (EEG) forward problems in anisotropic inhomogeneous media like the human head belongs to the class of the three-dimensional boundary value problems for elliptic equations with mixed derivatives. We introduce and explore the performance of several new promising numerical techniques, which seem to be more suitable for solving th...
متن کامل3D Thermoelastic Interactions in an Anisotropic Lastic Slab Due to Prescribed Surface Temparature
The present paper is devoted to the determination of displacement, stresses and temperature from three dimensional anisotropic half spaces due to presence of heat source. The normal mode analysis technique has been used to the basic equations of motion and generalized heat conduction equation proposed by Green-Naghdi model-II [1]. The resulting equation are written in the form of a vector –matr...
متن کاملA Lattice-Preserving Multigrid Method for Solving the Inhomogeneous Poisson Equations Used in Image Analysis
The inhomogeneous Poisson (Laplace) equation with internal Dirichlet boundary conditions has recently appeared in several applications ranging from image segmentation [1–3] to image colorization [4], digital photo matting [5, 6] and image filtering [7, 8]. In addition, the problem we address may also be considered as the generalized eigenvector problem associated with Normalized Cuts [9], the l...
متن کاملAdditive Operator Decomposition and Optimization-Based Reconnection with Applications
We develop an optimization-based approach for additive decomposition and reconnection of algebraic problems arising from discretization of partial differential equations (PDEs). Application to a scalar advection-diffusion PDE illustrates the new approach. In particular, we obtain a robust iterative solver for advection-dominated problems using standard multi-level solvers for the Poisson equation.
متن کاملA New Two-stage Iterative Method for Linear Systems and Its Application in Solving Poisson's Equation
In the current study we investigate the two-stage iterative method for solving linear systems. Our new results shows which splitting generates convergence fast in iterative methods. Finally, we solve the Poisson-Block tridiagonal matrix from Poisson's equation which arises in mechanical engineering and theoretical physics. Numerical computations are presented based on a particular linear system...
متن کامل